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Magnetically Accreting Isolated Old Neutron Stars 

Robert E. Rutledge 1 

ABSTRACT 

Previous work on the emission from isolated old neutron stars (IONS) ac- 
creting the inter-stellar medium (ISM) focussed on gravitational capture - Bondi 
accretion. We propose a new class of sources which accrete via magnetic inter- 
action with the ISM. While for the Bondi mechanism, the accretion rate Meondi 
decreases with increasing NS velocity, in magnetic accretors (MAGACs) Mmagac 
increases with increasing NS velocity (Meondi oc v~ 3 vs. Mmagac ^ v 1 / 3 ). MA- 
GACs will be produced among high velocity (^100 km s _1 ) high magnetic field 
(B >10 14 G) radio pulsars - the "magnetars" - after they have evolved first 
through magnetic dipole spin-down, followed by a "propeller" phase (when the 
object sheds angular momentum on a timescale ^10 10 yr). The properties of 
MAGACS may be summarized thus: dipole magnetic fields of -B^IO 14 G; min- 
imum velocities relative to the ISM of >25-100 km s _1 , depending on B, well 
below the median in the observed radio-pulsar population; spin-periods of >days 
to years; accretion luminosities of 10 28 -10 31 ergs -1 ; and effective temperatures 
£;T e ff=0.3 - 2.5 keV if they accrete onto the magnetic polar cap. We find no ex- 
amples of MAGACs among previously observed source classes (anomalous X-ray 
pulsars, soft-gamma-ray repeaters or known IONS). However, MAGACs may be 
more prevelant in flux-limited X-ray catalogs than their gravitationally accreting 
counterparts. 

Subject headings: pulsars: general - stars: magnetic fields - stars: neutron - 
X-rays: stars 
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1. Introduction 

Accretion onto isolated old neutron stars (IONS) through the spherical, Bondi-accretion 
mode ( |Bondi fc Hoyle 1944j ; [Bondi 1952| ) has been investigated as a means to produce a 



population of neutron stars detectable in the ROSAT X-ray All Sky Survey (RASS; [Treves 



|fc Colpi 19911 ; IBlaes fc Madau 19931 ; |Madau fc Blaes 1994j ; |Zaneet al. 19951 Popov et al. 20001 ) 



The estimated number of these sources in the RASS was initially high (~10 3 -10 4 ). However, 
as observations of higher mean velocities in the radio pulsar population ( [Lorimer et al. 1997 ; 



Hansen fc Phinney 1997] ; |Cordes fc Chernoff 1998| , and references therein) were considered 



the predicted number of detectable IONS accreting from the ISM decreased dramatically 
(~10 2 -10 3 ). This is because the Bondi accretion rate is a strong function of the NS velocity 
through the interstellar medium (ISM)Q - 



2 „„,-3 



M Bondi = 4vr(GM NS )> 

where p is the local mass density of the ISM. Bondi-accreting IONS are expected to have 
X-ray luminosities ~10 31 erg s -1 and effective temperatures fcT eff <; 200 eV, both dependent 
upon the accretion rate. 

There will be a different accretion mode, however, when the NS magnetic field dom- 
inates over the gravitational field in attracting accretion from the ISM. In this mode, the 
warm ionized medium (WIM) attaches to the NS magnetic field at the magnetosphere, ac- 
creting directly onto the NS surface. These magnetically accreting isolated old neutron stars 
(MAGACs, "magics") contrast with IONS which accrete through the Bondi-mode, in that 
the Bondi accretion rate decreases with the NS velocity, while in MAGACs, the accretion 
rate increases with increasing velocity: 



Mmagac = nR 2 u vp (2) 

where Rm = (p 2 /(4:irpv 2 )) 1 / 6 is the size-scale for the magnetosphere of a non-rotating NS, 
where the magnetic energy density from the NS's magnetic moment (p) is equal to the ram- 
pressure of the WIM. This gives Mmagac v 1 ' 3 . So, for a given magnetic field strength, 
Bondi accretion dominates at low velocities, while magnetic accretion dominates at high 
velocities. 



2 Here, v = \JV 2 + C^, the geometric sum of the spatial velocity of the NS and the sound speed in the 
accreted medium. The sound speed of the ISM is low (^10 km s _1 ) compared to spatial velocities relevant 
to the present work (v^lOO km s _1 ), so we assume ti>C s throughout. 
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The evolutionary scenario of magnetically accreting NSs has been outlined previously, 
by Harding & Leventhal ( |1992| , HL92 hereafter). HL92 found magnetised NSs first spin down 
via dipole radiation, followed by a longer spin-down period via the propellor effect. They 
restricted their attention to sources with surface magnetic fields ^10 13 G, and found that 
they are unable to accrete onto the compact object, as the propellor spin-down time is longer 
than the age of the universe. A similar line of argument is followed in the review by Treves 



et al. (|2000|) in the context of Bondi accreting IONS. This limitation does not apply to NSs 
with magnetic fields ^10 14 G; a NS with magnetic field B =10 15 G and velocity v =300 km 
s -1 , spins down in ~4xl0 9 yr. 

On the other hand, this same evolutionary scenario was used to argue for magnetic field 
decay in the IONS RX J0720.4-3125 QWang 1997| |Konenkov fc Popov 19971 ), and for NSs in 



general ( [Livio et al. 1998| ; |Golpi et al. 1998| ). We suggest, instead, that high magnetic field 



IONS have not been discovered due to the fact that strong magnetic field sources accreting 
from the ISM are spectrally harder (we find below) than Bondi accretors, and have been 
selected against by observers, or confused with background AGN. 

Moreover, since HL92's work, observations have revealed NSs with considerably stronger 
magnetic fields than considered by HL92, in the range of 10 15 G - the magnetars. Pulse- 
timing of soft-gamma-ray repeaters (|Kouveliotou et al. 1998| ; [Kouveliotou et al. 1999|) con- 



firms their existence as proposed in earlier theoretical work (Thompson fc Duncan 1993| ; 
Duncan fc Thompson 19921) . Circumstantial evidence suggests that anomalous X-ray pul- 



sars have similarly strong 5-fields ( Thompson fc Duncan 1996| ; [Vasisht fc Gotthelf 1997| : 



Heyl fc Hernquist 1997|) ; and precision pulsar-timing evidence also implies magnetic fields 



in excess of 10 14 G ( |Kaspi et al. 1999|) , confirming less-precise multi-epoch timing solutions 



which nonetheless led to a magnetar interpretation ( |Gotthelf et al. 1999| ). Static solutions 
for NSs have been proposed for surface fields up to ~10 18 G ( [Bocquet et al. 1995| ; |Cardall 



et al. 200T| ), although no objects with magnetic fields as high as this have yet been observed. 



Thus, with theoretical and observational motivation, we reconsider the scenario of HL92, 
but with particular attention paid to NSs with magnetic fields above 10 13 G. 

In §0 we describe the model for this population, following the scenario of HL92. In §H, 
we summarize the properties of MAGACs, discuss observational ramifications and conclude. 



2. Model 
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2.1. Evolutionary Model 

We assume that the neutron star is born with a velocity v = v 7 10 7 cm s -1 relative to 
the inter-stellar medium, and a magnetic moment /x = /J33 10 33 G cm 3 , which corresponds 
to a dipole field strength at the NS surface (10 km) of 10 15 G (appropriate for magnetars; 
[Thompson fc Duncan 1993| ; puncan fc Thompson 19921) . The spin period is initially short 



(P<C 10 sec). We assume the magnetic field does not decay. We assume a compact object 
mass M = Mns1-4 M , radius R = i?Nsl0 6 cm. We assume the NS travels through a spatially 
uniform proton density n p = 0.1n_i cm" 3 in the WIM (Poulares fc Cox 1990t [Ferriere 1998a| ; 
Ferriere 1998b| ). A more detailed examination of the properties of this population, which 



would include a distribution of NS kick velocities, orbits about a realistic galactic potential, 
and the galactic gas density and ionization state - such as that performed for Bondi accretors 
and cooling NSs ( [Popov et al. 2000| [Popov 200 1[ ) - will be the subject of forthcoming work. 



First, HL92 argues that magnetic accretion dominates over Bondi accretion when the 
magnetospheric radius Rm = 4.2 x 10 12 /i33 1//3 v 7 ^ 3 n_ 1 -1 / 6 cm is larger than the Bondi 
accretion radius -RBondi = 2G ^ NS ■ This inequality results in a velocity limitation for a 
MAG AC: 



> 94 /i 3 3- 1/5 n^ 10 M 3/ s 5 km s" 1 (3) 



Electromagnetic dipole radiation pressure inhibits accretion from the ISM until the light 
cylinder is external to the magnetospheric radius (-Rl.c. > -Rm), which occurs at a spin period 
of 

Pl.c. = 880 /i 33 1/3 v? 1/3 n_r 1/6 s (4) 
The timescale to reach this period through dipole radiation is: 

TL.c. = 2-3 x 10 7 J45 /^33" 4/3 vj 2/3 n_r 1/3 yr (5) 

where I 45 is the moment of inertia in units of 10 45 g cm 2 . This is a short timescale, compared 
with the life of the NS. After reaching this period, the light cylinder is at the Rm] however, 
the NS still cannot accrete (or, accretes inefficiently) due to the "propeller effect" ( [Illarionov 



|fc Sunyaev 1975|) in which matter is spun away by the magnetosphere, which is outside the 



Keplerian orbit in co-rotation with the neutron star (the co-rotation radius). This matter is 
ejected from the system, carrying angular momentum with it, spinning down the NS until 
the co-rotation radius reaches Rm, at which point the NS will have a spin period: 
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P P = 3.9 x 10 6 /i 3 3 1/2 vj 1/2 n_r 1/4 M NS 1/2 s (6) 

This spin period is substantially longer (by three orders of magnitude) than any spin 
period yet observed from a NS. The propeller spin-down timescale, when the accretion rate is 
dominated by magnetospheric accretion, is ([lllarionov fc Sunyaev 1975 ; Davies et al. 1979] ): 



r P = 2.3 x 10 10 J 45 /i 33 - 4/3 % 5/3 ™_r 1/3 yr (7) 

The Tp timescale dominates the evolutionary time between NS birth and emergence of the 
NS as a magnetic accretor. It is comparable to the age of the Galaxy (^10 10 yr) velocities 
already observed among the NS population, but only for high-magnetic field sources (^10 15 
G). This makes magnetic accretors with the right combination of magnetic field and velocity 
observable. 

Following spin-down via the propeller, matter is accreted onto onto the polar cap of the 
NS. The maximum accretion rate onto the compact object is (HL92): 

Mmagac = 7iR 2 M vp = 9.6 x 10 7 /i 33 2/3 v 1 ^ g s" 1 (8) 

This Mmagac is an upper-limit to the mass accretion rate onto the NS surface, and is 
potentially much lower (see § [2.2|) . Moreover, the Bondi and MAGAC accretion rates are 
simply related to the magnetospheric and Bondi radius: 

Mmagac ( Rm , , fV , 



of: 



Mpondi ^ -^Bondi - 

If we assume free-fall accretion onto the NS surface, this provides an accretion luminosity 



7 MAGAC 



1.8 x 10 28 /i 33 2/3 v)" n_! 2 / s M NS ergs" 1 . (10) 



To estimate the effective temperature, we take the size of the polar cap to be i? cap = 
Pns(#ns/Pm) 1/2 - This would seem unusually small - ~ 50 m for B =10 15 G. Nonethe- 
less, the effective temperature T e s of this emission is: 

kT eS = 0.39 /i 33 1/4 n_! 1/8 M^ s 4 R^ s keV. (11) 
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Note that kT e g is independent of the NS velocity. The accretion luminosity will remain below 
the local Eddington luminosity if: 

C"rs -^MAGAC f Rns\ 2 GMns ,~~\ 



tuh cttR 2 \ r J r 2 



where m# is the proton mass. The electron scattering cross-section in a strong magnetic field 
cx es is suppressed from the Thomson rate a? by <r es /<JT = 4 x 10 _6 (.E/1 keV) 2 (Bq E Y)/B) 2 , 
with 5qed=4.4x10 14 G, and photon energy E ( [Herold 1979| ). This result gives the constraint 



^33 -1 ^i 1/2 < 3 x 10 9 . This will not be violated for NSs with B >10 13 G unless the ISM is 
un-physically over-dense (n_i ^10 6 ). Thus, the accretion rate will be locally sub-Eddington. 

Finally, the ratio of Mmagac to Meondi is: 

MmAGAC = Lg ^2/3 ^10/3 n _ r l/3 M -| (1 g } 



Mi 



Bondi 



Thus, for the same velocity distribution, MAGACs will be more luminous than Bondi accre- 
tors. We discuss the observational ramifications of this in SB. 



2.2. Comments on the Accretion Luminosity 

We have two comments on the accretion luminosity. The first regards the accretion at 
magnetopause, and the second regards the need for the ISM to be ionized in order to be 
accreted. 



Arons & Lea ( 1976Q performed detailed analytic description of mass transfer across 



magnetopause in a 10 12 G NS system, fed by a Ri600 km s _1 wind from its companion, 
but with a substantially higher accretion rate at magnetopause (6xl0 19 g s _1 ) than we 
consider here. Even at this lower magnetic field strength and higher mass accretion rate, they 
find a cusp forms with a stand-off shock, causing an interchange instability which mediates 
accretion onto the compact object. While inflow at the cusp increases with increasing density, 
the interchange instability operates more efficiently and prevents very high densities from 



arising at magnetopause. Three-dimensional MHD simulations find similar results ([Toropin 
et al. 1999|). 



If the interchange instability of Arons and Lea ( |1976| ) permits accretion across the entire 
NS surface (rather than just at the magnetic poles), then the luminosity would remain the 
same, while dramatically decreasing the effective temperature, to: 
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kT cS = 0.0076 /i 33 1/3 v\ /l2 n^l & R^ 3 S /4 keV. (14) 

Further analysis of accretion onto a magnetic dipole with strength 10 13 -10 18 G is nec- 
essary to estimate the magnitude of the effect on the luminosity and spectrum. For now, we 
adopt Mm ag ac as an upper- limit. 

Secondly, to accrete material via magnetic interaction, the material must be ionized. 
Throughout, we have therefore assumed that MAGACs accrete from the (ionized) WIM, 
and not from the neutral ISM. However, one might wonder if the MAGAC's luminosity is 
sufficient to ionize the ISM by the time it reaches the magnetosphere (assuming that accretion 
is already taking place). Ionization of the ISM by an accreting IONS has been examined 
previously ( |Shvartsman 197T| |Blaes et al. 1995| ). These investigations found that a Bondi- 
accretion luminosity was capable of ionizing the ISM well beyond the Bondi-accretion radius. 
However, we have applied a similar approach to MAGACs accreting onto the magnetic polar 
cap, and find the fraction of a neutral ISM which would be ionized by the MAGAC at Ru 
to be <C 1, for a broad range of velocities (0.1< v 7 <10) and densities (0.1< n_i <100) for 
magnetic field strengths considered here. Therefore, MAGACs must accrete from the WIM, 
which has a typical density of ~0.13 cm -3 for a filling factor of 20% in the plane and a scale 
height of 1 kpc ( |Ferriere 1998b| ). This justifies our assumption that MAGACs must accrete 



from the WIM, and not from the neutral ISM. 



2.3. Limits on the Observability of MAGACs 



Here, we list those effects which will limit the observability of MAGACs, and which 
determine the luminosity, spectral hardness, velocity, magnetic field, and pulse-period pa- 
rameter space in which they may be observable. In Fig. |1], we divide the NS 5-field dipole 
strength, velocity parameter space, on the basis of these limits: 



1. MMAGAc>M"Bondi- As we find above, this results in a velocity limit (Eq. |3|), above 
which the NS is a MAGAC, and below which the NS is a Bondi accretor (broken line). 

2. The propeller phase must be shorter than the age of the galaxy. We use as a limit 
(Eq. 0) that t p <lxl0 10 yr (dotted line), and also include the line for r p <lxl0 9 yr. 
This age limit also effectively sets a magnetic field strength limit: NSs with B < 10 13 
G will not spin down in a Hubble time. Thus, if magnetic fields in NSs decay to <10 13 
G in less than a Hubble time, they will not be observed as MAGACs. 
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3. The pulsar must not escape the Galaxy. We adopt an escape speed limit of 700 km s _1 
(cf. |Miyamoto fc Nagai 1975] ; |Paczynski 1990| ), although the exact value depends on 
where in the galaxy the NS is formed. This limit is indicated by the hashed area at 
v >700 km s~ x . 

4. The NS magnetic field should not exceed the stability limit. We assume a surface dipole 
field of B <10 18 G (Bocquet et al. 1995| ; Cardall et al. 2001| ; dashed-dot line). 



We see in Fig. [j]that acceptable magnetic fields above 10 14 G, and acceptable velocities 
range from ~25 km s" 1 to the escape speed from the galaxy. We include lines of constant 
-^magac (for density = 0.1 proton cm -3 ; Eq. [10]); the luminosities of MAGACs are between 
10 28 and 10 31 ergs -1 (again, assuming free-fall accretion from the magnetosphere onto the 
NS). On the right axis of the figure, we mark kT e g, which is independent of velocity, and only 
weakly dependent upon the density of the ISM (oc n 1 / 8 ). MAGACs have /cT eff =0.3-2.5 keV. 



These are harder than observed from any of the identified IONSs to date (Treves et al. 2000| ), 
perhaps due to an observational bias. 

Finally, we note that for values of B =10 12 G and v =400 km s _1 which are not unusual 
in pulsars, such an object will not accrete via the Bondi-mode from the ISM, as the magnetic 
field is in the propeller range. In addition the source will not spin-down in <10 10 yr timescale; 
so it will wander the Galaxy as a dipole-spun down, propelling NS. 



3. Discussion and Conclusions 

Strong magnetic field (B >10 14 G) NSs go through three stages of evolution: (1) dipole 
emission; (2) propeller slowdown; and finally, the MAGAC stage (3) magnetic accretion onto 
the NS surface. MAGACs have the following properties: surface magnetic fields between 
10 14 G up through the NS stability limit (~10 18 G); spin periods between 10 6 -10 8 s; X-ray 
luminosities 10 28 -10 31 ergs" 1 , for a WIM density of 0.1 protons cm -3 ; and effective temper- 
atures /cT e ff =0.3-2. 5 keV. The low X-ray luminosities of these sources make it unlikely that 
the very long spin periods will be detected as pulsations with present X-ray instrumentation, 
as prohibitively long integrations with frequent monitoring would be required. 

At the median velocity for radio pulsars, MAGACs will be ~ xl30 more luminous 
than Bondi accretors (for v-j = 3, a typical median velocity for pulsars; [Lorimer et al. 19"97 ; 



Hansen fc Phinney 1997| ; |Cordes fc Chernoff 1998| , and references therein). Although the 



birth-rate of magnetars has been estimated to be ~10% of isolated pulsars flKouveliotou 



fc et al. 1994| ; [Kouveliotou et al. 1998| ), they can be observed to ~10 times the distance 
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and ~ 130 x the volume of the Galactic disk of Bondi accretors at the radio pulsar median 
velocity. Magnetic accretors, therefore, may be more prevalent in flux limited X-ray catalogs, 
such as the ROSAT All-Sky-Survey/Bright Source Catalog (RASS/BSC ; |Voges et al. 139^ ) 



than the lower-magnetic field Bondi accretors. Note, however, that Popov et al. ([Popov 



et al. 20001) found that the detected RASS/BSC INSs are more naturally interpreted as soft, 



young, nearby cooling neutron stars, rather than Bondi accretors; in this interpretation, no 
NSs accreting from the ISM have yet been discovered. The relative number of detectable 
MAGACs vs. Bondi accretors will ultimately depend on the particulars of the distribution 
of magnetic fields at birth in NSs, the kick velocities (and their possible dependence on B), 
as well as magnetic field decay, which we have neglected herein. 

A brief comparison between observations and the predicted properties of this population 
reveal no strong candidates for MAGACs among known source classes. 

Among radio pulsars, there are no known objects with implied B and v which would 
place them in the MAG AC B, v parameter space. The strongest B field radio pulsar is 
5.5 xlO 13 G ( |Camilo et al. 20"00| ), is just below the low-B limit; its velocity is unmeasured. 



At a velocity of ~ 300 km s -1 , this NS will spin-down on a timescale of rp ^210 11 yr, much 
longer than the age of the Galaxy. Thus, there are no observed radio pulsars in which the 
magnetic field is clearly strong enough to evolve into a MAGAC. This is possibly due to the 
relatively short timescale that such objects are dipole emitters (r L C ,^10 6 yr). 

The observed spin periods of magnetars are P ~1-15 sec, which, for the magnetic 
fields implied (10 15 G), are still within the range where electromagnetic dipole radiation 
will inhibit accretion. In addition, their X-ray luminosities are a factor ^1000 larger than 
those we expect for MAGACs. These objects therefore are not presently MAGACs; however, 
depending on their spatial velocities, they may evolve into MAGACs. 



The (model dependent) age of the isolated NS RX J185635-3754 (~1 My; |Walter 



et al. 1996| ; |Walter 2001]) is shorter than t p , unless the NS has close to B ~210 18 G. With 



such a high magnetic field, along with the implied velocity (?^200 km s _1 ), the source would 
be within the MAGAC parameter space. Its X-ray luminosity is ^2xl0 31 ergs -1 , which 
is within a factor of a few of the expected luminosity for a MAGAC with this B and v. 
However, the effective temperature (kT e s = 0.057 keV is substantially below that expected 
from a MAGAC with 5xl0 18 G (2.6 keV; Eq. [TTJ) . The spectral hardness seems to rule out 
this object as a MAGAC for the polar-cap accreting scenario. Since all previously observed 
INSs have effective temperatures < 200 eV ( [Treves et al. 2000|) , it seems similarly unlikely 
that any known INS is a MAGAC. It is possible, however, that an interchange instability 
at the magnetosphere permits accretion over a larger fraction of the NS, which results in a 
softer spectrum, such as observed here (ref. [L4|). 
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Discovering these objects in the RASS/BSC requires a similar approach taken in previ- 
ous work to identify IONS - identifying an X-ray source with no apparent optical counterpart, 
down to an X-ray-to-optical flux ratio of /x/Zopt^lOOO, with confirmation coming from an 
optical detection at an appropriate value. However, M AGACs can be spectrally harder than 
their Bondi accreting counterparts (kT e g ~ 1 — 3 keV; cf. [Popov 200 1| ), and searches for 



these objects must take this into account. Also, given their strong magnetic fields, these 
objects may have a cyclotron emission line, which may aid identification through X-ray 



spectroscopy ([Nelson et al. 19951) . Moreover, because we have assumed no magnetic field 
decay, the absence of MAGACs in the Galaxy, in combination with a population synthesis 
with high magnetic field SGRs and AXPs, would support field decay models. 

We conclude that MAGACs may present a population which are equally prevalent in 
flux limited X-ray catalogs as Bondi accreting sources. However, no observed examples of 
this class are yet known. 
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Fig. 1.— Magnetic field strength (B) vs. velocity for pulsars, showing regions where 
MAGAC IONS can be observed, and their accretion luminosities. The observability of 
MAGACs are limited at low velocities (left - dashed line) by a large magnetosphere (these 
then accrete always via Bondi); at high velocities (right - hashed area) by the likelihood they 
will escape the Galaxy's gravitational field; at low B (below - dotted line) from the fact that 
the NSs will not spin down via the propeller mechanism in <10 Gyr (we also include a line 
for Tp — 1 Gyr; and at high B (above - dashed-dot line) by the instability of high B neutron 
stars. Note also that a moderate magnetic field NS (10 12 G) with an average velocity (500 
km s" 1 ) will not accrete gravitationally due to the propeller mechanism. We include lines of 
the maximum accretion luminosity for the MAGACs after being spun down via the propeller. 
The right axis shows the effective temperature for accretion onto the magnetic cap, which is 
independent of velocity and only weakly dependent on density (oc n 1 / 8 ). 
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